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Due to the importance of fractional differential equations in solving many problems, researchers 
have paid attention to them in recent decades. One of the most important methods used to solve these 
equations is the hybrid spectral clustering method . 

We employed the law of impulse mass and Bernoulli's polynomials to solve this equation 


g (t, f(t), DP F(t), D*f(C), .., DP f (t)) =0,0<t<1fi)>fh,>-->>B, 
where the initial conditions are 

f"O)S%, t201; 2,105 =1, (2) 

In this research, we will rely on the fractional derivative in the sens e of Caputo. 

We will need Riemann-Liouville fractions to solve the fractional derivative and use numerical 


(1) 


methods to calculate it 


Preliminaries 
Now we will review the basic concepts that we will use in this research. 
Let D’ be the Caputo derivative operator and J’ be the Riemann -Liouville 
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fractional integral operator , fulfills the following property: 


m=—1 


[f-2pF FQ) —=D* FG) = D1 °O 7 ae 0<a<f,t>0, (3) 


ti-@ 
i=[a] 


The hybrid consisting of pulse mass functions and Bernoulli polynomials can be defined as 
follows: 


hin(t)n=1,2,..4.N  , m=T, 2,..,M, onthe interval /0, /) as 
n-1n 
B, (Nt - 1), te j—,= 
Inm(t) =} BmNe—n + 1) pal? 
0, otherwise, 


where vn is the order of bloc k - pu Ise func tions and mis the order of Bern oulli po lynomial Bm(t), 
which is given by the follo wing formula 


m 


Bit) = > @ Hocut™, 


k=0 
in which @ are the Bernoulli numbers. 


Riemann - Liouville integral of the hybrid functions 
The Riemann-Liouville equation is considered one of the most important hybrid equations 
using the fractional integration operator hnm(t), and the result was as follows: 


B n-1 n 
Drak ©)s N N 


B = is n 
Oe eG) be). 20S 


0, otherwise, 
where 


Lo are ae e @ mk or (1—nyk* 


n= rt+1 
r(r+1ythtr cP-1(" rt ¢) = Afnet 
T(r+B+1) (r+1)T(B) 2h (r+ 1L1— pir + 2; a N tr) ; 


and 
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k 


e0=> > 4 mene _ ayer 


k=0 r=0 


re tayerer tt (Fe) a tl” 
T(r+B+1) (+DI® iA (r+ 1a fir +2; 2(Fe)] 


where  F,(a,b;c;z) is the hyperg eometric function 
Suppose that 


H(t) = [hio(t), -.-» Rim (€), Rao (t), + Ram (0), Avo l(t), hvu (O)]"- 
Then, the Reimann-Liouville fractional integral of the basis vector H(t) is obtained as 
IPH(t) = HF (2), (4) 
whereas 


BOUL) =-[LP hag(t)yaiol? Rata Raglt pwegl” Hay (Piano gettin helt) 


Numerical method 
Now we will use the properties of the hybrid Bernoulli function that we mentioned previously to 
solve equations (1)-(2). To this end, let us assume that 


Do f(t) = ATH(t). (5) 


Where A is the unknown coefficient vector. using (3) with B = Bo and a = 0 and (3) , also taking 
the initial conditions (2) into account, we get: 


[Pol-1 i [Bol-1 ti 
f(t) =D F(t) + Y fOO)—= ATH) + DY fF (6) 
j=0 J j=0 J 
In the same way it works, fork = 1, 2,...,p, we get 
[Bol-1 : 
DP f(t) = AT HPo-Pe(t) + fi Pe. (7) 
’TG-Be +1) 


J=(Br] 
By substituting (5)&(7) into (1), we have 
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(Bolt, [Bol-1 i 
g| t, ATH%(t) + ». fi AT H(t), ATH0~#: (t) + >. fy ys ATF 0B () 
— j! r- [G=p. +) 

J=0 J=[B1] (8) 
[Bol-1 ti-Br 
+) fo J =. 
iG = 641 
a ee 
Suppose that 
es EU) goog 4 412.) FSO 
oy ee | ee ee 


which is the shifted Gauss-Chebyshev points. Now, consider t =¢,; in (8) a system of non 
linear algebraic equations in terms of the un known parameters of the vector A is obtained. 

We note that when solving the component system, we will obtain an approximate solution to 
problems (1) & (2) is given by (6). 

Example. We will now provide an example that demonstrates the accuracy and efficiency of the 
proposed method. We will choose the following fractional differential equations 


D°f() +f -VE-E=0,0<0<1 


f(0) =0 


The exact solution of this problem is f(t) = Vt . We have applied the method with 

different values of M and N. Table | Displays the absolute error at some specified points with M 
= 4 and N = 1, 4. Moreover, the numerical solutions obtained by N = 1 and M = 1& 2, together with the 
exact solution are plotted in Figure 1. 


Table 1. Numerical results with M = 4 and N = 1, 4 


r M=4,N=1 M=4,N= 
4 

0.2 5.982 x 105 6.231 x 
10° 

0.6 1.163 x 10% 9.682 x 
10° 

1.0 4.404 x 10° 3.674 x 
10-6 
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Figure 1. Numerical solutions obtained by N = 1 and M =1 & 2 together 
with the exact solution for Example 


Conclusion 

We have assumed a numerical algorithm based on the hybrid Bernoulli equation to solve High 
order differential equations. The algorithm relied on the fundamental properties of hybrid functions and 
transformed them into a nonlinear algebraic system. Then we solved an example using this algebraic 
system. 
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